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Abstract: We consider the problem of Arnold Diffusion for nearly integrable partially isochronous 
Hamiltonian systems with three time scales. By means of a careful shadowing analysis, based on a 
variational technique, we prove that, along special directions, Arnold diffusion takes place with fast 
(polynomial) speed, even though the "splitting determinant" is exponentially small. |^ 
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1 Introduction 

In a previous paper [| (see also §) we introduced, in the context of nearly integrable Hamiltonian 
systems, a functional analysis approach to the "splitting of separatrices" and to the "shadowing problem" . 
We applied our method to the problem of Arnold Diffusion, i.e. topological instability of action variables, 
for nearly integrable partially isochronous systems. The aim of this paper is to improve the shadowing 
theorem of Q and to apply this new theorem to the three time scale system ( |l.l[ ) below, in order to 
prove that along special directions Arnold diffusion takes place with "very fast speed" , namely a speed 
polynomial in e. To that effect, we use the results on the splitting provided in 

Three time scale Hamiltonian systems have been introduced in as a description of the D Alembert 
problem in Celestial Mechanics. Later on three time scale systems have been reconsidered for example 
in @, p, H, [0, §, |9|. 

In this paper we focus on isochronous three time scale systems as 



= -^/i + e"^/? • /2 + y + (cosg - 1)(1 + fifiip)), (l.I) 



where {ipi,(p2,q) G x T"^^ x are the angle variables, (/i,/2,p) € x R"~^ x are the action 
variables, (3 = {P2, ■ ■ ■ , Pn) € R"~^, n > 3, a > and e > 0, > are small real parameters. We will 
assume that /i = 0(min{e'^/^, e^°+^}). Hamiltonian describes a system of n isochronous harmonic 
oscillators with a Diophantine frequency vector uj^ = (l/-\/e, £°/3) , with one fast frequency uo^^i — \l\fe 
and (n — 1) slow frequencies lo^^2 = weakly coupled with a pendulum. 

When /i = the energy L0e,ili of each oscillator is a constant of the motion. The problem of Arnold 
diffusion in this context is whether, for fj, ^ 0, there exist motions whose net effect is to transfer 0(1)- 
energy from one oscillator to others in a certain time called the diffusion time. 

The existence of Arnold diffusion is usually proved following the mechanism proposed in |^. For 
/i = Hamiltonian admits a continuous family of n-dimensional partially hyperbolic invariant tori 
^0 = {(V'l 1:P) S T" X R" X X R^ I / = /q, q = p = 0} possessing stable and unstable manifolds 
W'{Ti„) = VF"(r/J ^ {{ip,I,q,p) e T"xR"xTixRi | / = /q, pV2 + (cosq-l) = 0} called "whiskers" 
by Arnold. For /z small enough the perturbed stable and unstable manifolds W^(7}^) and VF^(Tf^) may 
split and intersect transversally, giving rise to a chain of tori connected by heteroclinic orbits. By a 
shadowing type argument one can then prove the existence of an orbit such that the action variables / 
undergo a variation of 0(1) in a certain time called the diffusion time. In order to prove the existence 
of diffusion orbits following the previous mechanism one encounters two different problems: 1) Splitting 
of the whiskers; 2) Shadowing problem. 

^Supported by M.U.R.S.T. Variational Methods and Nonlinear Differential Equations. 



The "splitting of the whiskers" for Hamiltonian 7i^, when /i = O(e^), p > and £ ^ 0, has been 
studied in jl6|, and In and the size of the sphtting is measured by the "determi- 

nant of the sphtting matrix" which turns out to be exponentiahy small, precisely 0(exp(— (7r/2)e~^/^)). 
We underline that papers [p^-| 19 deal also with non-isochronous systems and more general perturbation 
terms (but two rotators only). 

In Q, the splitting of stable and uns table manifolds is related to the variations of the "homoclinic 
function" : T" — > R (defined in ( |2.5| )), which is the difference between the generating functions of 
stable and unstable manifolds at section {q — n}. VG'p(A) provides a measure of the distance between 
stable and unstable manifolds, so that a critical point A of gives rise to a homoclinic intersection. 
Usually det D^Gfi{A) is called the "splitting determinant". The use of the "homoclinic function" G^^ 
for measuring the splitting has two advantages. Firstly, it is very well suited to deal with the shadowing 
problem by means of variational techniques because G^ is nothing but the difference of the values of 
the Lagrangian action functional associated to the quasi-periodically forced pend ulum (2.2) at two true 
solutions, lying respectively on the stable and unstable manifolds VF^'"(7/q), see (2^). Secondly it may 
shed light on a "non uniform" splitting which would not be given by the splitting determinant, when the 
variations of G^ in different directions are of different orders. 

For the three time scale system associated to Hamiltonian 7i^, "non uniform" splitting is suggested 
by the behaviour of the first order expansion of G^ in /i, called the Poincare-Melnikov approximation. 



In fact the first order term, which is given by the Poincare-Melnikov primitive defined in (2.8), has 



exponentially small oscillations in the fast angle Ai, and polynomially small ones in the slow angles A2. 
Naively this hints the splitting to be exponentially small in the direction Ii and just polynomially small 
in the directions h- 

However, in general, for /x = 0(£^) and e — s- the homoclinic function G^ is not well approximated by 
the Poincare-Melnikov primitive. In iQ-l^ the asymptotic validity of Melnikov's integrals for computing 
the exponentially small "splitting determinant" is proved to hold only after exhibiting many cancellations. 

In § the naive Poincare-Melnikov approximation for Hamiltonian Ti.^ has been rigourously justified 
for iis~^^^ sufficiently small, in a different way. We define another "splitting function" G^, see (2.7), 
whose critical points as well give rise to homoclinic intersections. G^ is well approximated, for fj, = 0(£^) 
and e — !■ 0, by the Poincare-Melnikov primitive and has exponentially small oscillations in Ai, see theorem 



2.2. The crucial observation is that Gu and Gu are the same function up to a diffeomorphism %fjn of the 



torus close to identity, namely G^ = G^ o 7/;^, see theorem 2.1. 

After the works [|),§,|2§l,|2),|l§,§,[|l3),@and references therein, it is a well estabhshed fact 
that the diffusion time is estimated by a polynomial inverse power of the splitting. For instance, using 
the estimate on the size of the splitting of and an exponentially long diffusion time has been 
obtained in namely Td — 0(exp(G/e'')) for some 6 > (see also theorem 5.2 of |^). 

However the properties of G^ (oscillations of different amplitude orders according to the direction) 
suggest that Arnold diffusion can take place with different speed along different directions; since, for 
larger splitting one would expect a faster speed of diffusion, one could guess the existence of diffusion 
orbits that drift along the "fast" directions I2 G R"^^, where the splitting is just polynomially small 
w.r.t. 1/e, in a polynomially long diffusion time = 0(l/£''). The aim of th is paper is to prove that 
this is indeed the case. In order to prove this phenomenon (see theorem for the general case and 
theorem for an application) we refine the shadowing theorem 2.3 of |^ for dealing with the present 
"non-uniform" splitting. Note that, because of the preservation of the energy along the orbits, Arnold 
diffusion can take place in the direction I2 for 7i > 3 only. 

In order to justify heuristically our result we recall how the diffusion time Td is estimated in ^J, 
once it is verified that stable and unstable manifolds split. Td is, roughly, estimated by the product of 
the number of heteroclinic transitions k (= number of tori forming the transition chain = heteroclinic 
jump/splitting) and of the time Tg required for a single transition, namely Td = kTg. The time for a 
single transition T^ is bounded by the maximum time between the "ergodization time" Tf, of the torus 
T" run by the linear flow w^t, and the time needed to "shadow" homoclinic orbits for the corresponding 
quasi-periodically forced pendulum equation 2.2. 

The reasons for which we are able to move in polynomial time w.r.t 1/s along the fast I2 directions 



are the following three ones, (i) As in M, since the homoclinic orbit decays exponentially fast to 0, 



the time needed to "shadow" homoclinic orbits for the quasi-periodically forced pendulum (2.2) is only 
polynomial, (ii) Since the splitting is polynomially small in the directions I2, we can choose just a 
polynomially large number of tori forming the transition chain k = 0(l/eP) to get a 0(l)-drift of l2- {Hi) 
Finally, the most difficult task is getting a polynomial estimate for the "ergodization time" -defined as 
the time needed for the flow {io^t} to make an a-net of the torus- with a appropriately small. By a result 
of Q this time satisfies Tg = 0{l/a'^). Let us explain how this estimate enters into play. In order to 
apply our "gluing" variational technique, the projection of our shadowing orbit on the torus T", namely 
{tUct + Aq}, must approach, at each transition, sufficiently close to the homoclinic point A to be capable 
to "see" the homoclinic critical point A of G^. The crucial improvement of the shadowing theorem 4.1 



allows the shadowing orbit to approach A only up to a polynomially small distance a = 0{s^), p > 0, 
(and not exponentially small as it would be required when applying the shadowing theorem of |^). By 
the forementioned estimate on the ergodization time Te — 0(l/a^) it results that the minimum time 
after which the homoclinic trajectory can "jump" to another torus is only polynomially long w.r.t 1/e. 
Actually this allows to improve as well the exponential estimate on the diffusion time required to move 



also in the Ii direction, see remark 4.3 



Theorems [4.1] and |4.2| are the first step to prove the existence of this phenomenon also for more general 
systems (with non isochronous terms and more general perturbations). 

The paper is organized as follows: in section 2 we recall some preliminary results taken from |^]. 
In section 3 we introduce the general "splitting condition" which will be used in section 4 to prove the 
shadowing theorems. 

Through the paper Ci and 6i will denote positive constants which are independent of s and fi. 



2 Preliminaries 

In this section we recall the results of that will be used in the sequel. We refer to Q for complete 
details and for the description of the general functional analysis approach based on a Lyapunov Schmidt 
type reduction. With respect to the notations of Q we remark that we have changed the sign of the 
perturbation / in Hamiltonian 7i^. 

The equations of motion derived by Hamiltonian are 

(p^uje, i = -n{cosq-l) d^f{(f), q=P, p = sinq {1 + ^ f{(p)). (2.1) 



The angles (p evolve as ip{t) — uj^t + A; therefore equations (2.1) can be reduced to the quasi-periodically 
forced pendulum equation 

-q + smq{l + ^lf{uJ^t + A)) = 0, (2.2) 

corresponding to the Lagrangian 

Cf^^Aiq, g, t) = y + (1 - COS q)il + fifiiu.t + A)). (2.3) 



For each solution q(t) of (2.2) one recovers the dynamics of the actions I{t) by quadratures in (2.1). 

For /i = equation ( ^.2[ ) possesses the one parameter family of homoclinic solutions to 0, mod 27r, 
qe{t) = 4 arctan(exp (t — 9)), 6* G R. Using the Implicit Function Theorem one can prove (lemma 2.1 
of 1^) that there exist, near the unperturbed homoclinic solutions qe{t), for < /i < /xq small enough 
independently of ui^, "pseudo-homoclinic solutions" g{t) of equation ( |2.2| ). These are true solutions of 
( |2.2D in each interval (—00, 9) and {9, -l-oo); at time t = 9 such pseudo-solutions are glued with continuity 
at value g(9) = n and for t ±00 are asymptotic to the equilibrium mod 27r. We can then define 
the function i^^ : T" x R — > R as the action functional of Lagrangian (|]^) evaluated on the "1-bump 
pseudo-homoclinic solutions" g{t), namely 

C^AqA eit), ql e t)dt+ C^Aq'X e € eit), t) dt, (2.4) 
-oo Je 



and the "homoclinic function" : T" — > R as 

G^{A) ^ F^{A,0). (2.5) 

There holds 

F^{A,9)^G4A + u,d),ye en. (2.6) 

Remark 2.1 The homoclinic function G^ is the difference between the generating functions j^{A, q) 
of the stable and the unstable manifolds VF^'"(7}u) (which in this case are exact Lagrangian manifolds) 
at the fixed section {q — tt}, namely G^(v4) — S~ j^{A,n) — S'^ j^{A,n). A critical point of Gf^ gives rise 
to a homoclinic orbit to torus 7}^, see lemma 2.3 of j^. 

In order to justify the dominance of the Poincare-Melnikov function when /i — O(e^) one would need 
to extend analytically the function F^(A, 6*) for complex values of the variables. Since the condition 

g{Re 9) = tt, appearing naturally when trying to extend the definition of g to € C, breaks 
analyticity, the function _F)j(A,0) can not be easily analytically extended in a sufhciently wide complex 
strip. To overcome this problem, in |^ the Lagrangian action functional is evaluated on different "1- 
bump pseudo- homoclinic solutions" g. Define -00 : R- ^ R- by '!/'o(*) = cosh^(t)/(l + cosht)'^ and set 
il^git) = ipQ{t — 0). Two important properties of the function '!/'o(0 are that /p^ V'o(^)'?o(^) dt ^ Q and 
that it can be extended to a holomorphic function on R + i(— 7r,7r) (while the homoclinic solution q^it) 
can be extended to a holomorphic function only up to R + i(— 7r/2, 7r/2)). By the Contraction Mapping 
Theorem there exist (lemma 4.1 of [|j) near qe, for /z small enough, pseudo-homoclinic solutions Q'^ g{t) 
and a constant g defined by 

-Qi^^g + sin Qi^/l + f{cJet + A)) ^ ai^^gMt) and j^{Q%{t) - qe{t))Mt) dt ^ 0. 



We define the function : T" x R ^ R as the action functional of Lagrangian (2.3) evaluated on the 
"1-bump pseudo-homoclinic solutions" Q^g{t), namely 



F^{A,e)^ f C^,A{Q^^git),Q'Xg{t),t) dt 

Jr 



(2.7) 



and Gf, : T" ^ R as Gf,{A) = F^(A, 0). 



Remark 2.2 Also critical points of Gfj^ give rise to homoclinic solutions to torus Tj^, see lemma 4-2 of 
By theorem below, from a geometrical point of view the introduction of the "homoclinic function" 
G^ may be interpreted simply as measuring the splitting with a non constant Poincare section, see the 
introduction of 

The crucial point is now to observe that the homoclinic functions G^ and are the same up to a 
change of variables close to the identity, as stated by the following theorem (see theorem 4.1 of [^) 

Theorem 2.1 For fi small enough (independently of uJe) there exists a Lipschitz homeomorphism (a 
real analytic difjeomorphism if f is analytic) "tjjn : T" T" of the form -tpulA) = A + kfi{A)uj^ with 
fc^ : T" ^ R satisfying fc^(A) = 0{fi), |fc^(A) - fc^(A')| = 0{n)\A - A'\ such that = o ^p^. 

Let r(e, A) denote the Poincare-Melnikov primitive 

r{e,A)= [ {1 - cos qo{t))f{u,t + A) dt. (2.8) 



/R 

Develop in Fourier series w.r.t. the first variable the homoclinic function G^{A) = X^fciez 5^1 (^2)6*'"'^ '^^ 
and the Poincare-Melnikov primitive r(e. A) = Ffc^ (e, A2)e''^i'^^ . Assume that the perturbation / 

is analytic w.r.t {(p2, ■ ■ ■ , fn)- More precisely assume that there exist > for i — 2, . . . ,n, such that / 
has a C°° extension in _D := R x (R-t-«[— r2, r2]) x . . . x (R + i[— r„]), holomorphic w.r.t. {(p2, . ■ ■ , ifin)- 
Denote the supremum of |/| over D as ||/|| := sup^g^, |/(<<5)|. The following theorem about the splitting 
of stable and unstable manifolds in three time scale systems, holds (see theorem 5.1 of ||]) 



Theorem 2.2 For small there holds, for all Ai £ T"--^ 

G^(Ai,^2) = goiA2) + 2Re[gi{A2)e'^']+RiAi,A2) 

= Const + (fjro{e,A2) + Ro{e,ti,A2)) +2 Rc \fiTi{e, A2) + Ri{e, fi, A2) 



RiAuA2), 



where Ro{e, /i, A2) ~ 0[ii 



Ri{e,n,A2) = O 
RiAi,A2) = o[^^e-^/^ 



■ exp 



and 



exp 



In order to prove our shadowing theorem we need also to recaU the definition of the k-bump pseudo- 
homoclinic solutions q\ g (t) for the quasi-periodically forced pendulum ( 2^ ) . Such pseudo solutions 
turn k times along the separatrices and are asymptotic to the equilibrium 0, mod 2tt, for t — > ±00. More 
precisely in lemma 2.4 of |^ it is proved that for all fc G N, for all 6*1 < . . . < 0^ with min^ (6*^+1 — 9i) > L, 
with L sufficiently large, independently of LUg and ^, there exists a unique pseudo-homoclinic solution 

g{t) : R ^ R which is a true solution of ( |2.2| ) in each interval (—00, di), {Oi, Oi+i) (« = 1, . . . , fc — 1), 
{9k, +00) andqj^gie,) = 7r(2z-l),g^g(i) = q'^^g^it) in {-^Ji) and q^git) = 2TTk + q^g^{t) in (Ok, +00). 
Such pseudo-homoclinic orbits arc found via the Contraction Mapping Theorem, as small perturbations 
of a chain of "l-bump homoclinic solutions" q^g - 

Then we consider the Lagrangian action functional evaluated on these pseudo-homoclinic orbits q^ g 
depending on n -|- A; variables 



-t-00 



^^MAe{t)AAfi{t),t) dt. 



Setting efe — (!,...,!) G R'^, the following invariance property, inherited from the autonomy of Ti^, holds 



Let T^, : T" X R'^ 



F^{A, 9 + Tjek) = Fl:{A + tjuj,, 9), V0 G R^ e R. 
R be the "fc-bump heteroclinic function" defined by 
J^'^{A,9) ■.^Fl:{A,9)-iIl,-Io)-A. 



(2.9) 



(2.10) 



Lemma 2.1 V/0,/0 ^ H^"", if {A, 9) is a critical point of the "k-bump heteroclinic function" T^{A,9), 
then {If^{t),u!^t + A, qA,g{t),q^g{t)) where ly^if) = Iq - ^ /^^(cosq^ g(s) - l)d^f{uJeS + A)ds is an 
heteroclinic solution connecting Tf„ to Tji^ . 

By lemma |2.l| , in order to get heteroclinic solutions connecting 7/„ to 7}/ , we need to find critical 
points ofT^^{A,9). When mmi{9i+i — 9i) +00 the "fc-bump homoclinic function" F^'(A, 6*) turns out 
to be well approximated simply by the sum of the functions F^{A, 9i) according to the following lemma. 

We set 9o = —00 and 9k+i = +00. 

Lemma 2.2 There exist positive constants Ci,Li > and functions Ri{ii, A,9i_i,9i,9i^i) such that 
Ve > 0, VO < ^ < ^0, VL > Li, V6'i < . . . < 6*^. with mmi{9^+i - 9i) > L 

k k 

F^^{A, 01, • • • , 0,0 = E 0,) + E R^{^i, A, 9,.i,9,, 0,+i), (2.11) 

i=l i=l 



with 



\R,{fi, A, 9,^1,9,, 0,+i)| < Ci exp(-Cii). 



(2.12) 



3 The splitting condition 



We now give a general "splitting condition" on the homoclinic function well suited to describe the 
non-uniform splitting of stable an d u nstable manifolds which takes place in three time scale systems. 
Roughly, the "splitting condition" 3.1 below states that Gf^ possesses a maximum and provides explicit 
estimate s of the non-uniform splitting. It will be used, in the next sec tion, to prove the shadowing 
theorem |4.l|. As a paradigmatic example, we will verify, in lemma 3.2, that, when the perturbation 



f{f) — X^'i^i COB if j, the "splitting condition" is satisfied, see also remark 3J 



Condition 3.1 "Splitting Condition". There exist A G R" and a basis {ili, . . . ,^n} o/R", n > 3, 
such that oj^ G R+fJi, 1/2 < |r2j| < 2, detjrii, . . . , r2„} > 1/2, {fis, . . . , f2„} is an orthonormal basis 
of {f^i, ^2}'^ , and which enjoy the following properties : let us define iJ^(ai, . . . , a„) as the homoclinic 
function G^{A) in the new basis, namely 



Gf^iA + aifli 



(3.1) 



Then there exist positive constants p, cr, (5i, (52, <^3 > 0, with 3cr < p, 62 < 5^, and two continuous functions 

2 n — 2 

hj2 '■ [—PiP] X B — > R with li{x) < l2{x) for all x G [—p,p\ x B , such that: 



(i) /or X = (02, ... , a„) G [-p, p] x B^ 



■n-2 



J{x):= sup H^{ai,x)>TCL&yi\H^,{li{x),x),H^^{l2{x),x)\+8i\ 

aie[h(x)M(x)] ^ ^ 



(a) for all y = (03, 



,an)eB , 



Va2G[-a,a], J{a2,y) > J{0,y) - ^, 
Va2 G -p + 2a] U [p - 2a, p] , ^(02, y) < ^(0, y) ~ 62; 



(Hi) 



7i — 2 

Va2 G [-a, a], \/y eB^ , 



J{a2,v)>J[().Q)-^. 



Va2G[-p,p], VyGS; \B;zI, J{a2,y) < J{0,0) ~ S,. 

The next lemma states that the former "splitting condition" is satisfied by the homoclinic function 
G^ if (and only if) it holds for the homoclinic function G^. 



Lemma 3.1 Assume that Gf_i satisfies the splitting condition 3.1 with maps li,l2 and parameters p, a, 
Si, S2, S3. Then G^i satisfies the splitting condition 3.1 as well, for some maps li^2 = li,2 + 0{fi/ \/e) and 
with the same parameters. The converse is also true. 

Proof. By theorem G^^ = G^ o ip^, where ipfj.iA) = A + ^^(A)^^^ and i/'/^ is a homeomorphism. Set 
Hf_,{ai, ...,an)^ Gf_i,{A + ai^li -|- . . . -f- a„17„). We have 



Hf^{ai,a2, . . . ,an) = Hf,[ai + fc^(ai, . . . ,a„)— — -,02, 



where k^{ai 



, an) := kf,{A + aiQi + . . . + a„r2„)- 



Assume that satisfies condition 3.1 with maps h^h- For all x — (02, 



G [-P, p] X Bp , 



the map ai 1-^ ai + kf^{ai, x)\lo^\/\^i\ is a homeomorphism from the interval (/i(a;), ^2(2;)) to the interval 



{Ii{x),l2{x)), where lj{x) :— lj{x) + kfj.{lj{x),x)\uJe\/\i^i\ {j = 1,2). There results that, for all x 



n — 2 

(a2, . . . ,a„) e X B 



J{x) -.^ JVl^^ Hf^{ai,x) ^ sup Hfj,{ai,x) ^ J{x). 



Therefore satisfies the splitting condition 3.1, with maps Ij replaced by Ij, and the same positive 
parameters. Since fc^ — 0(/i) and = 0{l/y/e) we have \lj{x) — lj{x)\ = 0{^i/y/e). ■ 

We now give a paradigmatic example where the former "splitting condition" is satisfied. Assume 
that the perturbation / is given by f{^pi, . . . , (fin) = Sj=i costpj. In the next lemma we show that the 
corresponding homoclinic function satisfies the "splitting condition" 3.1 and hence, by lemma 3.1, G^ 



as well satisfies the "splitting condition" 3.1 



Lemma 3.2 Assume that f{ip) — X]j=i cos ipj . There exist a basis {^i, . . . , ri„} and a positive constant 

Sq such that, if e is small, < iie^^^^ < Sq and < /is^^"^^ < Sq, then satisfies the "splitting 
condition" IJ, with A = 0, p = Tre^+i/^ a = p/6. Si = S3 = pp^/2, 82 = Snpe"^/^ e^p{~Tr/{2y/e)), 
li{x) = -2ii,l2{x) = 2tt. 

Proof. In order to simplify the notations we give the proof for ri = 3 and we assume that — 1. We 
will prove that satisfies the "splitting condition" 3.1 with ^ = and w.r.t the basis 



r!i = (l,e"+i/2/5), 



(0,/3), 



^3 



(0,/?'), 



where |/3'| = 1 and /3 • /?' = 0. We set p — Tre^+^Z^ and we assume that < /i < Sp^ , < /i < fe'^ /^, where 
(5 is a small constant (independent of e) to be specified later. Let 5 > be such that theorem 2.2 holds 
for < /X < Se^/"^. We shall always choose Q < S <S. 

From now on, notation Ki will be used for positive universal constants, whereas notation Ci{S) will 
be used for positive constants depending only on 5. Notation u = 0{v) will mean that there exists a 
universal constant K such that lul < K\v\. 



Our first aim is to prove expression (3.14) below. It easily results that, if /((p) — X]j=i cos ^Pj, 

3 



ro(e, A2) = .A„i,il^a^\ cos and ri(e, A2) 



U sinh(/3,£-f ) 



By thereom 2.2 we have 



G^ (^1,^2, A3) =30(^2,^3) + 2Re 51(^2,^3)6^^^ 



and, by (|3.2D, up to a constant that wc shall omit, 



30(^2,^3) = 



_J^;^^M^cosA2 
sinh(/32e'^f ) ' sinh(/33e ^ 



51(^2,^3) 



2 > 



-7r/2v^ 



Vesinh(5^) 
In this proof we shall use the abbreviations 

27r/32£° 27r/33£° 



)■ 



(3.2) 

(3.3) 

(3.4) 
(3.5) 



G, 



sinh(/32£"f) sinh(/33e°f)' 



27r 



Note that, as e ^ 0, we have 
27r/3j£'' 



/esinh(^) 

• ^=4 + 0(£"), D.^^e-^/'^^^'Ul + Oie-^ 

smh(/3j£'»|) V£ ^ 



(3.6) 



We shall consider e small so that 



(3.7) 



By (|1|) and (|1|), since < < 

G,, ( Ai , yl2 , A3 ) = go (^2 , ^3 ) + M^e cos Ai + O f 



-l/2p-7r/v^_ 



-1/2-^/2^ 



(3.8) 



Since {Ai, A2, A3) — oiSli + a2^2 + 03^^,3 = (oi, (aie^+^Z^ + 02)/? + a3/3'), the homoclinic function in 
the new basis {f2i, $72, ^^3} writes 



Hf,{ai,a2, as) = G^ (aif^i + a2f72 + a3f^3) = (ai, (aie^+^Z^ + 02)/? + 03/?') . 



(3.9) 



Define /io(fo2, 03) = 50(^2/? + a3/3')- By ( |3.9| ) and (3.8), there exists cq{S) > such that, for all < e < 
co{S), 

H^iai,a2, as) = ho (aie^+^Z^ _^ 03) + cos ai + o(/ifc"i/2g-^/2v^j ^3 -^q^ 



We derive from this latter expression and (3/7) that 



H, 



,(ai, 02, 03) = /lo (aie°+'/' + 02, 03) + o(^^ie-'^^e-^/^^ 



By (^ and (|3J) 

50(^2,^3) = A^G, - 2f,{Al + Al) + 0(iie'^(Al + Al)) + o(/i(A4 ^ ^4^^ ^ o(^2)_ 



(3.11) 



(3.12) 



We shall assume in the sequel of the proof that 02, 03 G [— p, p], ai € [— 27r, 27r], so that, since p = 7re°+^/^, 
there resuhs ai£''+^/^ € [-2p, 2p], 62 = aie"+^/^ + 02 S [-3p, 3/o] and 63 + a| = 0(/). Moreover we have 
that < ^(5p^ and there exists Ci{5) £ (0,co(i5)) such that, if < £ < ci{S), then e"" < 6 and p'* < Sp"^. 
Note also that, since (3, (3' are orthonormal vectors, we have A^ + A\ = + a^. Finally we derive from 
(|3l|) that, for < e < ci((5). 



ho{b2, as) ^ fiC, - 2p.{bl + al) + 0{p.5p^). 



(3.13) 



Since p = Tre""*"^/^ we have e ^/^e '"1'^^ — o{p^) as e — + 0; therefore, by ( 3.11 ) and ( |3.13| ), there exist 
Kq > 0, C2((5) e (0, ci((5)) such that, for all < e < €2(6), 



Hfj,{ai,a2,a3) = pCe ~ 2p{h\ + 03) + ro(ai, a2, 03), |?'o(ai, a2, 03)] < Kofj,Sp^ 



(3.14) 



where 62 = oie^+^Z^ + 02. 

We now prove that point (i) of the "splitting condition" 3.1 is satisfied by G^ with 61 = p,p^ /2, li{x) = 
— 27r and l2{x) = 27r where x := (02, 03). Let us consider J{a2, as) '■— supjj^gj_27r 27r] H^{ai, 02, 03). Since 
02 G [— p, p], — a2e^'-"+"'^/^'' G [— TT, tt] and we can derive from ( 3.14 ) that 

J{a2.as)>H^(^-a2e-^''+^'^\a2,as) > pC, ~ 2iial - R'opSp^ (3.15) 

If ai = ±2tt then 62 = ai£"+-'^/^ + 02 = 02 ± 2p and then, since 02 G [— p, p], we get I&2I > P- As a 
consequence, by ( p^ and (|1|), 



( ± 27r, 02, 03) < pGe - 2/i(p2 + al) + KopSp^ < J(a2, 03) - 2pp'^ + 2Kop5p^ 



(3.16) 



Cho osing 5 < l/2Ko, we get in ( 3.16 ) that H^{±2Tr, a2, a3) < 1/(02, as) — pp /2. It results that condition 
3J-(i) is satisfied with Si = pp^/2, /i(a;) = — 27r and ^2(2;) = 27r where x = (02, 03). 



We now turn to the proof of (ii) and (iii). If 02, 03 £ [—p,p],ai £ [— 27r,27r] and|62| = \aie°-~^^^^+a2\ > 
V2K^p, then, by ( p^ and ( p5| ), 



Hf^{ai,a2, as) < fiC^ ~ 2/1(03 + 2A'o(5/9^) + K^piSp^ < pC^ - 2pa\ - ZKopSp^ < ^7(02, 03). 



Hence 



^^(02, as) = sup {-ff^lai' as) ; aie°'+'^''^ G - 02 - ^J'i.KoSp, -02 + \JlK^bp | 



(3.17) 



We use here that, since 2Kq8 < 1, [-02 - V^K^p, -aj + V2i?;)7p] C [-2/?, 2p] = [-27re''+i/2^ 27re'^+i/2] 
Writing ai = (62 - a2)£"("+^/2\ we derive from (|3.17D and (|3.10D that 



sup 



b2e[-\/2Ko5p,\/2i?o5p] 



sup 



HA{b2 - a2)e' 



02,03 



0(62, a3)+M^e cos (^H_^))+0(<5/ie-i/2e--/2v^). (3.18) 



b2 6[-v'2Ko5p,\/2i?o5p] 

Now, if 62 e [-\/2^o^p, V2^o^/o] then 626"""^/^ G [-7rV2^o?, 7r\/2^o^], so we can write that 

'^"^^1 =cos (^^^ 



cos 



(|^)=»(?375)+«(^)- 



(3.19) 



As a consequence, by ( 3.18 ) and (3/7) there holds 



^^(02, 03) 



where we have set 



sup 



m 



[h,{b2, a,) + pD, cos (^^)) + 0[V~6ps-'/'e~-/'^) 



(03) + pD, cos (^^) + 0(^VSpe-' 



-l/2g-7r/2v^ 



771(03):= sup ho{b2,a3). 

62 G [- \/2Ko^p, \/2Ko5p] 



(3.20) 
(3.21) 



FinaUy, there exists Ki > such that, by ( 3.2Cl| ), 
:7(a2,a3) =^(0,a3)+A*^e(cos(-^^)-l)+ri(a2,a3), |ri(a2,a3)| < K^VSpe-^/^e-^^^^ . (3.22) 



We are now in position to prove condition 3.1-(m). Assume < S < ir'^/AKf and choose a = p/6 = 
ga+i/2^/g and 62 = inpe-^/^e-''/^^. IfaaG [-a, a] then cos(a2/e''+^/2) - 1 > -1 + \/3/2 > -1/6. This 
readily implies, by (|3.22[) and (js^), that ^'(aa, a3) > ^^(0, 03) - S2/2. If 02 G [-/?, -p + 2a] U [p - 2g, p ] 
then 02/6"+^/^ e [-TT, -27r/3] U [27r/3,7r], so that cos(a2/e°+^/^) - 1 < -3/2. It follows, stiU by (|3.22| ) 

and (p/7|), that 1/(02,03) < J^( 0, 3) — (52. This proves cond ition 3^-(m). 

In order to prove condition 3A_-{iii), we notice that, by ( [3.13| ) and the definition of m given in ( ^.21 ), 



"^(03) — pCe — 2pa\ + 0{5pp^). 



(3.23) 



Hence there exist K2 > Q and C3(J) G (0,C2(^)) such that, for aU < £ < C3((5), by ( |3.23| ), ( |3.20| ) and 
(13) 

^^(02,03) =^(0,0)- 2/102 + ^2(02^^3), with 1^2(02,03)1 <i^2W- (3.24) 



Let us assume 5 < I/6K2 and let S3 — pp^/2. By ( 3.24 ) and (3/7), if 03 e , then ^'(02,03) > 

J(0, 0) - <53/2; if 03 e bI'\b';zI, then J(o2, 03) < /7(0, 0) - 5^. 

As a conclusion, lemma 3^ holds with (5o = min{(5, 1/2Kq^'k'^ /AKf, l/6i^2}. ■ 

Remark 3.1 The former splitting condition holds also also for more general perturbations f{fi, . . . , fn) 
for which /o(</52, • • ■ , ^n) possesses a nondegenerate maximum at (^2j • ■ • j^n) '^'^^ /i(^2; ■ • ■ ;^n) 7^ 
where /fc^ (932, • . ■ , V^n) = (l/27r) /q /(u, </32, . . . , </5n)e~*'"'^'^ dcr. T/iis kind of condition is considered in 
theorem 5.2 of m. 



4 The shadowing theorem 



In this section we shall prove, under the "splitting condition" 3.1, our general shadowing theorem. 



Theorem 4.1 Let n > 3 and assume that the homoclinic function satisfies the splitting condition 3.1 



Let Lo^ be a (7^, T)-diophantine vector, i.e. \lu^ ■ k\ > 'y^/\k\'^ V/c g Z"\{0}. Then, for all /q, Lq £ R" such 
that (/g — /q) € Spanjris, . . . there exists an heteroclinic trajectory from Tj^ to Tj'^ which connects 

a rj -neighbourhood of torus Tj^ to a rj-neighbourhood of torus Tj'^ in the "diffusion time" 

Td<C^ti^ max|^,|ln5i|,|ln52|,^| + |lnry|, (4.1) 
w/iere A := max , , — n-2 ^2(2;) — min , , —n-2li{x)\. 

Remark 4.1 The diophantine condition on the frequency vector lo^ restricts the values of e and (5 that 
we consider. In any case, if for instance f3 is (7,71 — 2) -diophantine then for r > n — 1 there exist Cq > 
and a sequence Sj — > such that ui^ is (7^, t)- diophantine with 75 = co£°, see for example 



Remark 4.2 The meaning of (^.1) is the following: the diffusion time Td is estimated by the product of 
the number of heteroclinic transitions k = ( heteroclinic jump / splitting ) = — /ol/f^a, and of the time 
Ts required for a single transition, that is Td ~ k ■ Tg. The time for a single transition Ts is bounded by 
the maximum time between the "ergodization time" i.e. the time needed for the flow cut to make 

an a-net of the torus, and the time max{| In^ij, | ln<52|j A/lwgj} needed to "shadow" homoclinic orbits for 
the forced pendulum equation. We use here that these homoclinic orbits are exponentially asymptotic to 
the equilibrium. 

We could prove also the existence of connecting orbits for all 1q — Lq G Span j^h, ■ • ■ , ^n}- In this 
case the number k of heteroclinic transitions would depend also on 62, see remark ^.S. 

Proof. Still for simplicity of notation we write the proof for n = 3. Then /q — /q = ±|/q — /ol^a; we 
assume for definitiveness that Iq ~ Iq — \Iq — /o|f^3, so that (/q — /o) • {J2'j=i "j^j) = l^o ^ -^o|a3- 
We choose the number of heteroclinic transitions as 



r8|/^-/o|p 



S3 



1. (4.2) 



By lemma 2.1, in order to prove the theorem, it is sufficient to find a critical point of the fc-bimrp 
heteroclinic function : T^ x R'^ R such that 

e,-e,= o( ^'^° ~^°' max{ I ln<5i|, | ln52|, ^|) ■ (4.3) 

We introduce suitable coordinates (oi, 02, a^, si, . . . , Sk) G R'^ x {—2tt, 2tt)'' defined by 

A = A + J2'^j^i and Vz = l,...,fc, = i_!]l±Il_^^})\^ ^ (4.4) 
j=i 1'^"' 

where rji are constants to be chosen later. Let H^{a, s) = FI^{A, ff) be the "A;- bump homoclinic function" 
and H^{a, s) — ^fHA, 9) be the "fc-bump heterochnic function" expressed in the new variables (a, s). 
The function Ti.^ does not depend on ai, since, by the invariance property (2^) (we recall that 

= |17iK/|c^e|), 

KM = F^,{A + ±a,n,, fa + , . . . , + '^;7^^'"^' ) - (^0 - /o) ■ {ta,n, 

= F^, {A + ±a,n,, fa , . . . , fa "^' ) - K - Io\a3. (4.5) 



In the sequel of the proof we shaU use the abbreviation TYj^ ~ ^^(0, 02, 03, s)- 

We now choose the constants (771, ... , 77^) e R''. Note that, since cj^ is (7^, T)-diophantine, fii satisfies 
the diophantinc condition 

|f^,..|>^, v.ez"\{0}. 

Hence, by the results of [Q, there exists C > such that the "ergodization time" Te of the torus T'^ 
run by the linear flow fJii, i.e the smallest time for which {ilit ; < t < Te\ is a a— net of the torus, 
can be bounded from above by C\ijJ^\I {^^a'^). Hence for each interval J of length greater or equal to 
C|a;e|/(7e(T'^) there exists 77 € J such that 

d{Tflx,2-KT?)<a. (4.6) 
In particular there exists a constant C2 and there exist 7/j such that 

7^ In {^C,—^-—\ + A < 7,+, - ^7. < ^ m (SC, . + ^ + A, (4.7) 

6i|i2i| V mm{di,(32}/ Ci|S2i| V mm{di,(32}/ 7etT^ 

77iili EE Xi, mod27rZ^, Xi = 2/4^2 + Zi^a with |yi| < cr, |zi| < cr. (4.8) 



In order to prove the theorem we just need to prove the existence of a critical point of TC^ in x 
(minZi,maxZ2)'^- The upperbound of the diffusion time given in (4.1) will then be a consequence of ( [4.7D 
and {^). Indeed, by (^^ and (^) we get that 



7i+l — Ci 



iVt+i - Vi)\^i\ , (si+i - Si)l^il ^ J_ 



,ln(. 



8C1- 



1 



{61,62}^ 7eCr^ 



(4.9) 



By (4.9) there exists C > such that the time 6*^+1 — 9i "spent for a single transition" is bounded by 

Ts := max(0,+i - 0,) < Cmaxl | ln(5i|, | ln<52|, (4.10) 

I ^"feCr Pel J 



From (4.10) and ( [4.2| ) we derive immediately (4.3) and then ( [4.l[ ). 

We now provide, using lemma 2.2, a suitable expression of the fc-bump heteroclinic function TC^. By 
lemma 2.2, the invariance property (|2.6D, (4.8) and since : ^ R, we get 



H^(a2,a3,s) = ^ [F^(A + ^ajl7j, 

i=l j=2 

k _ 3 



{r], + s,)\ni\ 



OJs 



I Qt \ 1-^0-^0 1 

+ bi[a2, as, Si-i,Si, Si+ij ; 



El' 

i=l 
k 

n' 



J=2 



-03 



\i!^-io\ 



as + Si 



(4.11) 



where S'i := S'i(a2, as, Si_i, s^, s^+i) = i?i(yl., 6'j, 6'j+i) after the change of variables (4.4). The left 
hand side inequality in (4.7) implies that 



Ji+l — Pi 



> 



im+i - r^, - A)\ni\ ^ 1 



> — In 8C1 — 



hence, by (2.12), 



\S,\ < 



Ci 

min{(5i, (52} 



min{(5i, 62} 



(4.12) 



We will maximize in the open set 
[/ = |(a2,a3,s) e R''+^ Vz 02+?;^ G (-p, p), as+z^ G (-p, p), e (li{a2+yi,a3+Zi), 12(0-2+^1, a3+Zi)jY 

U ^ ^ since {0} x {0} x Il^^i{li{yi, Zi), hivi, Zi)) C U. Since U is bounded, attains its maximum over 
U at some point (a, s) = (02,03, s). It is enough to prove that (a, s) € C/. 

• We first prove that for alH, Si e ('i(a2 + yjjtis + ^i); ^2(^12 + 2/11 ^3 + -s^i))- Since (a,s) is a maximum 
point of in U, for any t e [li{a2 + 2/^,03 + ^i), ^2(^2 + Vii^s + Zi)], replacing Si with t does not 
increase Ti.^. Since such a substitution alters at most three terms among Si, . . . , Sk in (4.11), we 
obtain, using (4.12), that for any i, for any t G [h{a2 + yi,a3 + Zi), ^2(^2 + ^1,^3 + Zi)], 

3 

Hf^{s.i,a2 + yt,a3 + Zi) > -ff^(t, 02 + y^, 03 + Zi) - - min{Si,S2}. 



Hence 



,_ _ _ x _ X 3(5i 
J^f,(si, a2 + yi, a3 + Zi) > J (02 + yi,a3 + Zi) — 



and, by condition 3. l-(i), this implies that G (Zi(a+Xj), Z2(a+Xi)), where we have set Xi ~ {yiTZi). 
We now prove that for all i, 0,2 + yi G (— p, p). Indeed we have by ( 4.1l| ) and (4.12) 

H^(a2,a3,s) < [g.(^. «2 + aa + ^0 + '""^^^^^ - ^^^3 



< 



(4.13) 



On the other hand, still by ( 4.11 ) and ( 4.12 ), choosing s = s = (si, . . . ,Sfc) with Si G (^1(2/1, 03 
Zi),l2iyi,a3 + Zj)) so that H^{si,y^,a3 + z.,) = J{yi,a3 + z^), we get 



:(0,a3,S) = ^[i?. 



m1 5ijyii«3 + Zi 



\I',-h\ 



03 + 5"^ 



i=l 
k 



> 



> 



2^ [J{yi,a3 + z^)- — as 



i=l 
k 



(4.14) 



since \yi\ < <J and by condition |3.1[-(m). Since H^ia,!) > H{^(0,a3,s), we can derive from ( 4.13 ) 
and ( 4.14 ) that there exists iq G {1, . . . ,k} such that ^7(02 + yio,a3 + Zig) > J'{0,a3 + Zi„) — 3(52/4. 
Still by condition ^]^(m), 02 + yia G (— p + 2cr, p — 2a). As a result, since (by (4^)) \yi — yio \ < 2cr 
for all i, we get that a2 + J/i G (— p, p). 



• At last we prove that, for all i, 03 + z^ G (— p, p). By (4.11) and (4.12), choosing s — s = (si, . . . , Sk) 
such that Si G {li{yi, Zi),l2{yi, Zi)) so that H^{si,yi, z^) = J'(yj, z^). 



H^(0,0,?) = E[i/, 



i{(5i,<52} 



^E['^(0>0)-^], (4.15) 



since ^2 < S3, \yi\, \zi\ < a and by condition 3.1-(zm). Hence, since H^{a,s) > 7^^(0,0,3), by ( 4.13| ) 
there exists zq G {1, . . . ,k} such that 



J{a2 + y.„,a3 + z,,) + I - ^^^S3 > J(0, 0) - ^. 



(4.16) 



Since jasl < p and by we get |/^ - /oUaal/A: < 5^/8. Hence by ( |4.16| ), J'(a2 + j/i,,, as + ZiJ > 
J7(0, 0) — 753/8. By condition 3.1-(«m), 03 + G (— p + 2<7,p— 2a) and as a consequence, since for 
all i < cr, we deduce that 03 + Zi S {—p,p). 

We have proved that the maximum point (a, s) S J7, which completes the proof of the theorem. ■ 



As a consequence of the general shadowing theorem 4.1 and of lemma 3.2 we get the following theorem 



Theorem 4.2 Let f{ip) — fj! " ^ 3, and cj^ 6e a (7^, T)-diophantine vector. Assume e, p.e 



-3/2 



and /ie 



to be sufficiently small. Then, for all with oje ■ lo = We • I'q and (/o)i = (^o)i there 



exists a heteroclinic orbit connecting the invariant tori Tj^ and Tii with a diffusion time 



Td<C- 



ps' 



7777:77 X max , 

a+(l/2) ^Jei£ 



(4.17) 



Remark 4.3 The number k of heteroclinic transitions used in the proof of theorem J^.i is polynomial 
w.r.t 1/e since our shadowing orbit moves along the directions {I2, ■ . ■ , In) G R"^^ only ("directions of 
large splitting"). On the contrary the shadowing orbit connecting toriTj^ andTj'^ with (/q)i 7^ (-^0)1 would 
shadow an exponentially large number of heteroclinic transitions and the diffusion time would depend also 
on the contant 82 — 37r/ie~"'^/^ exp(— 7r/(2y^)) which describes the exponentially small splitting. In any 
case, at each transition, the shadowing orbit approaches the homoclinic point only up p = 0(e"+^/^) and 
therefore the time Tg spent for each single transition is polynomial w.r.t 1/e. In this way we deduce that 
the diffusion time Td is estimated, up to inverse powers of 1/e, by an exponential Td = 0(exp(7r/(2Y^))). 
Since the (determinant of the) splitting A — 0(exp(— 71/(2^^))) we get that Td ~ while in JT^ and 

^ the diffusion time is estimated by Td ~ l/A^ for some positive constant p. 
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